General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



c^-nj oij 


AgRISTARS 


available under NASA sponsorship, ^ 
In the interest of early and wide dis- 
semination of Earth Resources Survey 
Program information and without liability 


for any use made thereoi." 

Supporting Research 




SR-Xl -04033 
NAS9-15981 

E83'10309 

A Joint Program for 
Agriculture and 
Resources Inventory 
Surveys Through 
Aerospace 
Remote Sensing 

January 1981 



ASSEMBLY LANGUAGE CODING FOR CLASSY 

(EB3-10309) ASSEMBLY LANGUAGE CODING (fUR N83-27299 

CLASSY (Elogxc, Iiic.) 36 p HC AU3/Mf' AJl 

CSJL U2C 

Unclas 
G3/43 00309 


M. E. Rassbach 



Elogic, Inc. 

4242 S.W. Freeway, Suite 304 
Houston, Texas 77027 



Lyndon B. Johnson Space Center 

Houston. Texas 77058 


1, Report No, 2, Government Accession No. 

SR-X^ -04033 

3. Recipient's Catalog No. 

4. Title and Subtitle 

Assembly Language Coding for CLASSY 

5. Report Date 

6. Performing Organization Code 

7, Author(s) 

11. E. Rassbach 

8. PBiiorming Organization Report No, 

NAS811 

10, Work Unit No. 

9. Performing Organization Name and Address 

Elogic, Inc. 

4242 S.W. Freeway, Suite 304 
Houston, TX 77027 

11. Contract or Grant No, 

NAS9-1 5981 

13. Type of Report and Period Covered 

Technical Report 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Lyndon B, Johnson Space Center 
Houston, TX 77058 

14. Sponsoring Agency Code 


IS. Supp4tm«ntary Notes 


16. Abnract 


A set of assembly language improvements to the CLASSY clustering 
algorithm have been developed and tested. These are descrived 
in detail and analyzed. 


17. Key Words (Suggested by Author(sM 18. Distribution Statement 

CLASSY algorithm 

19. Security Qassif. (of this report) 20. Security Classif. (of this page) 21. No. of Pages 22. Price* 

Unclassified Unclassified 35 


JSC Form 1424 <R«f Nov 75) 


For sale by the National Technical Information Service, Springfield, Virginia 22161 


NASA — JSC 




















SR-X1-O4033 
NAS 9 -15 9 81 


ASSEMBLY LANGUAGE CODING FOR CLASSY 

BY 

M. E. RASSBACH 


This report describes Classification activities of the 
Supporting Research project of the AgRISTARS program. 


Elogic, inc. 

4242 S.W. Freeway r Suite 304 
Houston, Texas, 77027 


Jeuiuary 20, 1981 



CONTENTS 

Purpose 1 

Matrix Arithmetic 2 

New Exponential Routine 5 

Path Difference in CLASSY 9 

Appendix 1: Program listings for vector routines . 11 

Appendix 2: Program listings for the fast 

exponential system 26 


Assembly language coding for CLASSY 


Elogic has developed several assembly language routines to 
speed up the CLASSY algorithm. 

These fall into two catagories 

1) matrix arithmetic 

2) improvement in the exponential routine 


Purpose of assembly-coded Matrix Arithmetic and Exponential 
Routines . 

Currently available compilers produce a great number of 
extra instructions/ particularly "overhead” instructions such 
as Load/ Store, and index arithmetic. These instructions can 
be eliminated by careful hand coding of the routines. Thus it 
is possible to reduce the machine-processing time for a program 
by coding the most frequently used routines in assembly language. 

CLASSY is particularly suited for this form of optiaization, 
since most of its execution time is spent in a few vector and 
matrix processing routines. In addition, as was noted while de- 
signing the vector processing routines, application of the distri- 
butive law will allow putting one multiplication outside of the 
inner loop in several instances. (It was later realized that this 
change could be made in the Fortran version, as well.) Thus by 
using new, assembly language versions of the matrix arithmetic 
routines, CLASSY could be greatly accelerated. , Some results are 
given in Table 1. 
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Matrix Arithmetic 

The matrix arithmetic routines correspond exactly (bit-by- 
bit) to Fortran versions ("new Fortran”) of the old matrix 
routines. The new Fortran versions differ from the original 
Fortran versions ("old Fortran") by the application of the 
distributive law and precalculation of some quantities to reduce 
the number of operations inside the inner loop. Rounding 
error differences can make the new version differ from the old 
version, which may become substantial due to the characteristics 
of the CLASSY algorithm (See "Path Differences in CLASSY" below) . 
However, since these differences are only due to floating-point 
rounding error, which is random in both the new and old versions, 
the results are equivalent. The new version may in fact give 
slightly higher precision, due to the method of summation used. 
Table 2 gives the routine names for the new and old Fortran 
versions, as well as the assembly version. The names of EXEC'S 
used to switch from one version to cmother are also given. The 
names given are the neunes of the source Icmguage files; the ob- 
ject progrcuns must have the original name. 
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Time 


Version 

Routines 
(Vector, XP) 

Iterations , 
Channels 

Time 

(Improve- 

ment) 

per Cluster 
per Iteration 
(Improvement) 

Original 

Old Fort., Fort. XP 

10 

4 

369 


4.10 


Old Fort. 

Old Fort., ALC XP 

10 

4 

518 

(0.71) 

3.99 

(1.03) 

New Port. 

New Fort. , ALC XP 

10 

4 

301 

(1.23) 

3.34 

(1.23) 

ALC 

ALC, ALC XP 

10 

4 

240 

(1.54) 

2.67 

(1.54) 

Original 

Old Port., Fort XP 

7 

8 

1590 


15.14 


Old Fort. 

Old Fort., ALC XP 

7 

8 

1164 

(1.37) 

18.47 

(0.82) 

New Fort. 

New Port. , ALC XP 

7 

8 

927 

(1.72) 

9.46 

(1.60) 

ALC 

ALC, ALC XP 

7 

8 

641 

(2.48) 

6.54 

(2.31) 


Table 1: Run time comparison. 

This table shows a comparison of the run times of CLASSY 
using different combinations of routines. The run was 10 iterations 
on 4 channels of the simulated data, without any map generation. 
Initialization and randomization time has been subtracted off. Im- 
provement is the speed ratio relative to the corresponding original 
version. (ALC means Assembly Language Coding, Port, means Fortram.) 

All runs were done on a simulated data set using a 4900 
point data set. The scatter in the results is due to the effects of 
path differences (see heading on these) . 
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Table 2: Routine ncunes for matrix arithmetic. 


Routine 

Old Fortran 

New Fortran 

Assembly 

DOTSQ 

DOTSQB 

DOTSC 

DOTSQA 

CORECT 

CORECTB 

CORECT 

CORECTA 

MPVS 

1 

MPVS3 

MPVS 

MPVSA 

VPV 

VPV 

WV 

VPVA 

VfCTV 

VMTV 

VMTV 

VMTVA 

EXEC to compile 

VERSOLD 

VEl^NEW 

VERSALC 
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New Exponential routine (XP) 

In order to speed up executions, especially when the 
nuntber of cheumels is few, Elogic has written a special ex- 
ponential routine^ called XP, specialized to the needs of the 
CLASSY algorithm. In CLASSY, the exponential is calculated 
as often as is the primary quadratic form, and the standard 
exponential routine requires several multiplies and a divide. 
The new exponential used for CLASSY reduces these times by 
table look up and by allowing a small increase in the RMS 
error. The new routine lias an RMS precision of about 1 part 
in 24,000, which is easily adequate for the needs of CLASSY. 

In addition, the new routine uses linear interpolation and 
requires only one multiply. 

The new exponential routine calculates XP (x) » exp (-x/2) 
by subdividing the interval between 1 and 16 into 240 equal 
parts. A table contains the coefficients of the least-squares 
fit of a linear form (First - order (polynomial) to the desired 
exponential in this region. These coefficients give a result 
accurate to about one part in 24,000 RMS in the result. Argu- 
ments outside the range of 1 to 16 are handled by directly ob- 
taining the exponential of the remaining part and multiplying 
it times the part calculated above. Note that because of limi- 
tations on the exponent reuige of the result to the routine, 

176 

arguments are limited to lie between 2 log 10 . This routine 

will not raise an exponent underflow condition if the result is 
too small; rather, it returns 
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a result of zero. This saves additional time since STATIS 
will not be required to check the arg\Mnent range before cal- 
culating the exponential. 

The table used XP is currently calculated at execution 
time by the routine XPREP which must be called by CLINIT 
during initialization. (A version of CLINIT making this call 
must be used whenever the XP routine is to be used. ) XP could 
be modified to contain these constants directly. A modified 
version must be also used of STATIS , which calls the XP 
routine. 

The new exponential routine XP(x) = exp (-x/2) is thus 
a fast assembly Icuiguage routine designed to calculate the 
oft-repeated exponential function quickly and with an accuracy 
easily sufficient for use in CLASSY. 

The XP routine causes slight variations in the calculated 
probabilities in CLASSY, which can cause the algorithm to taka 
a different path (See the note "Path Differences in CLASSY" 
below. ) 



ORIGINAL PAGE IS 
OF POOR QUALITY 
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XP routine 

The exponential e~^^^ is calculated as follows : x is divided 

into three parts, I, n, and E. 

X » 161 + ^ + S 

-22 < I < 21 ; 0 < n < 256 ; 0 < 2 < 5 — 


Then 

2 XPBIG(I)*XPAR(n) * (XPD + 2) 
z 0 I < -22 or I > 21 

This approximation has a relative RMS error of about 1/27000, 
which is well adequate for ithe CLASSY statistical system. The 
variablCJ OUTCNT counts the number of overflows made. 

XPP, XPAR, and XPBIG are generated by the routine XPREP. 

Equations for constants used in exponential routine: 

XPAR(n) = d^e"^“”"“ n = 0(1)256 

XPD = 

d- = sinh a/a c 1 + + • • • 

^ ® 2 

di = (sinh a - a cosh a) ~ J ^ . 

81 * 1 “X /2 

a “ T “ ^ intervals of in original argiiment (of e ' ) 

These coefficients are the result of a least-squares fit of 
(dg + d^ 2 ) to for ^ < 2 < I 

XPBIG (I) = e”®^ “23 £ I < 22 

These are arranged in common block XPCOM (which is part of 
the exponential routine) as follows: 
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Quantity 

Type 

Bytes 

OUTCNT 

1*4 

4 

XPXTRA 

1*4 

4 

XPD 

R * 8 

8 

XPAR 

R * 4 

1024 

XPBIG 

R * 4 

354 


XP procedure; 

The argument mantissa is shifted into the lower half of a 
double word (CARVE ) , by adding a special constant. If the argument 
X is 0 £ X < 16, the upper word is fixed, and is checked. The 
upper byte (n) of the man .i.., ia is used as an index to XPAR, and is 
then replaced with a standardr sign-exponent byte, to yield C. (The 
non-normalized character of 5 does not matter because it is immed- 
iately added to XPD, which has a larger exponent. This also elim- 
inates the random 4 low-order bits of 5.) The calculation is then 
direct. For x outside (0,16), the procedure is the same, except 
for a range check and an additional factor of XPBIG, derived from 

an index calculated from the upper half of the double word. If 

+76 

X is out of range (+352, corresponding about to 10— ) , the result 

is set to 0, and overflows are counted. 

CARVE structure; 


mantissa of Z . 



decimal point 


I 
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Path Differences in CLASSY 

The CLASSY algorithm may take widely different paths 
to the maximum likelihood solution, depending on very small 
effects. In some cases, whenever the clusters have not been 
fully resolved statistically due either to too few iterations 
or too little data, the two paths may end at slightly different 
points. Variations causing the algorithm to go on different 
paths include rounding error (for example use the associative 
and/or distributive lav^s) use of new standard function routines 
(such as the exponential routine), or otlier minor changes. 

The reason CLASSY is so susceptible to small changes in 
the current version is the use of ^iOnta-Carlo methods in pro- 
cessing the points with respect to classes which give them low 
probi&biiity. Small fluctuations in the calculations may accumu- 
late until the threshold for the Monte-Carlo process is exceeded 
in one version and not in another. Since the random number in 
generator is sequential, this offsets the usage of all the suc- 
ceeding ramdom numbers, giving completely different Monte-Carlo 
choice for each point. 

This, in turn, leads to moderately large fluctuations in 
the numerical values used, which ultimately changes the point at 
which the SPLIT or other threshold is passed. Since one cluster 
being SPLIT and separated cam make a major difference in the pro- 
cessing of other clusters, the path taUcen by CLASSY rapidly 
becomes quite different. 
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In summary, minor changes in CLASSY can make large differences 
in the path the algorithm takes. Tlxis is mediated by the amplifi- 
cation caused by various thresholds, in particular those used in 
the Monte-Carlo subsystem. All the paths should, however, 
converge to the same point, except when insufficient iterations 


or data are used. 
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Appendix 1: Program listings for vector routines 

This appendix contains listings of all the routines 
listed in Table 2. They are in the order: old Fortran, 

new Fortran, and assembly language. Preceding the routine 
listings is a listing of the EXEC routines for converting 
one set to the other. 



ORIGINAL PAGE IS 
OF POOR QUALITY 
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SUBROUTINE CORECT <REL,PV,P,S) 

CONNON /RISC/ NQ,NH,LR,LU,NINCL3,HXAR,UTINIT,KROOT,EPS,OELT, 

1 AHa,ODCON,XOUFLO,XUNFLO,li'AIIJIN,ELll1TH,SEPTH,VFAC,AHH,SBLTH, 

2 INDXVL,UFAC,NPTSO,PQRATH,SPHUTH,DUFAC,GRACTN,ANOFAC, 

3 AH0HIN,AH0HAX,AN0RAT,V0LLIH,BIAS,PJ0IN,VRJ01N,USIN,UDELSH, 

4 BETTER, NODE, C0RLEN,SPC0R 
REAL REL(30), PV(30), S(30) 

DO to 1 > 1,H0 

REL<I) » PV<I) - S(I) / P 
II » I 

C URITE (4,9999) II,REL(I),PV(I),S(I),P 
9999 FORHAT ("CORECT I,REL,PV,S,P",I4,4(E10.4,2X)) 

10 CONTINUE 
RETURN 
END 


COR00010 

COR00020 

COR00030 

C0RO0C40 

COROOOSO 

COR00040 

COR00070 

COROOOSO 

COR00090 

COR00100 

C0RO0110 

C0R00120 

COR00130 

COROOHO 

COROOtSO 
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ORIGINAL PAGE IS 
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SUBROUTINE CORECT (REL,PV,P,S} 

COHHON /HISC/ NQ,NH,LR,LV,NINCLS,HXAR,UTINIT,KROOT,EPS,DELT, 
t AHQ,00C0N,X0UFL0,XUNFL0,UADJIN,ELIHTH,SEPTH,UFAC,ANf1,SBLTH, 

2 INDXUL,UFAC,NPTS0,PQRATH,SPHVTH,DgFAC,6RACTH,AH0FAC, 

3 AHOHIN,ANONAX,AHORAT,VOLLIN,BIAS,PJOIN,URJOIN,USIH,UDELSH, 

4 BETTER, NODE, C0RLEN,SPC0R 
REAL RELI30), PV(30), S(30) 

C 

P0V*1./P 
DO to I > t,NQ 
REUI) = PV(I) - S(I)*P0V 
C Uf^ITE (A,99?9) I,REL(I),PV(I),S(I),P 
9999 FORNAT (■'CORECT I,REL,PV,S,P'',I4,4(EI0.4,2X)) 

10 CONTINUE 
RETURN 
END 


COROOOlO 

COR00020 

C0R00030 

COR00040 

COR00050 

C0RO0P60 

COR00070 

coRoooao 

COR00090 

COR00100 

C0RO0110 

C0RO0120 

COR00130 

C0RO0140 

C0RO01SO 

C0R00160 
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CORECTA 

•* THIS ROUTINE UAS WRITTEN BY ELOGIC,INC., N0V,I?7?, TO REPLACE THE COR00010 
* CORRESPONDING FORTRAN ROUTINE IN THE CLASSY SYSTEM. FOR ADDITIONAL COR00020 


:* INFORMATION CONTACT ELOGIC 

,INC. 

COR00030 



COR00040 

!> SUBRGITINE CORECT(REL,PV,P, 

S): REL(I)*PV<I)-S(I)yp 

COR00050 

CORECT CSECT 


COR00060 

USING CORECT, 15 


COR00070 

STM 2,?, SAVE 


COR0OO8O 

LM 1,4, 0(1) 

ARCS ADDR 

COR00090 

LA 8, 4(0,0) 

INCREMENT IS 1 

COR00100 

EXTRN MISC 


COR00110 

L 7,=A(MISC) 


COR00120 

USING MISCD,7 


CQR00130 

L 9,MQ 

COMPARAND IS MQ 

COR00140 

SLA 9,2 

MULTIPLY BY 4 

COR00150 

SR 9,8 

SUBTRACT 4 

COR00160 

■4 


COR00170 

LE 6,=c'\.d' 


C0R00180 

OE 6,0(3) 

GET 1/P 

COR00190 



CQR00200 

SR 5,5 

DO LOOP INDEX I 

COR00210 

IC LE 2, 0(4, 5) 

S(I) 

C0R0O220 

MER 2,6 

MULTIPLY BY 1/P 

CQR00230 

LE 0,0(2, 5) 

PV(I) 

COR00240 

SER 0,2 

SUBTRACT REG2 

COR00250 

STE 0,0(1, 5) 

STORE REL(I) 

CQR00260 

BXLE 5,8, IC 


CQR00270 

:» 


COR00280 

LM 2, 9, SAVE 


COR00290 

BR 14 


COR00300 

SAVE DS 9F 


C0R00310 

4 COMMON /RISC/ 


COR00320 

MISCD DSECT 


COR00330 

HQ OS F 


CGR00340 

END 


COR00350 


of PO°" 


WlGINAL PAGE IS 
OF POOR QUALITY 
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FUNCTION D0T5Q(V,AHET) 

CALCULATES THE INNER PRODUCT U.V RELATIVE TO THE METRIC ANET 

DINENSION L1ST<38), XTHP( 8), YTHP( 7), UADJ(24), U(26), 

t INDEX(37), LSUBS<40),LSUPER(39), 1DADJ(23), NSYNB<36>, 0U(25), 

2 PCUN<33),PRIRCN<32), CIN(22), CT0T(21), PRQP<2?), SPFAC(tS), 

3 0PR0P(23), V0LIN<18), V0LRT(17), DC0N<U), PQRAT(14), 0DEN(19>, 

4 0ISS(35), PPASS<34), PST(30), 0CIN(2O), PC0ND(3t ),0PRI0R(27), 

5 PAVE<t3), PILE(12) 

DINENSION VRIN(475),GEN<999),GREF<999),ALINK(99) 

EQUIVALENCE (LINK(41 ) ,ALINK(4I ) ,GREF(8) , GEN(7) ,VRIN(7) ) 
EQUIVALENCE ( LINK(41), LSUBS<40),LSUPER<39), LIST<38), 


INDEX<37), NSYHB(36)» 
PRIRCN(32), PC0ND(31), 
OPRIOR(27), U(2A), 
CIN<22), CT0T(2t), 


DISS<35), PPASS(34), PCUN(33), 

PST(30), PR0P(29), OPROP(28), 

3 OPRIOR(27), U(2A), 0U<25), UADJ(24), I0A0J(23), 

4 CIN<22), CT0T(2t), 0CIN<20), 0DEN(19), VQL1N(18), 

5 V0LRT(17), DC0Nn4), SPFACttS), PQRAT(14), PAVE<t3), 

6 PIL£(12), XTNP( 8), YTHP< ?)) 

COHNON/CLUS/ JUNK< 1 2) ,NARL,NT0P,NTBSZN,NUANT,LINK(14000) 
DINENSION NXAR(31),LR<3),LV(3) 

EQUIVALENCE (LR<1 ),LVRIN),<LR(2),LKURT), 

\ <LR(3),L0VAR)r(LV<1),LSUN),(LV<l),LSKEU),<LV(3),L0SUN) 

CONNON /NISC/ NQ,NN,LR,LV,NINCLS,NXAR,UTINIT,KROQT,EPS,DELT, 

1 ANQ,ODCON,XOVFLO,XUNFLO,UADJIN,ELINTH,SEPTHtVFAC,ANN,SBLTH, 

2 INDXVL,UFAC,NPTSO,PQRATH,SPNVTH,DUFAC,GRACTN,ANOFAC, 

3 AN0N1N,AN0NAX,AN0RAT,V0LLIN,BIAS,PJ0IN,VRJ0IN,USIH,UDELSN, 

4 BETTER, NODE 

< END OF STANDARD BLOCK. ) 

DINENSION NTB<32) 

CONNON /STPAR/UAIT,C0NLV,SKBND,SKCHI,TRBN0,TRCHI,URK6ND,URKCHI, 
1 PACCEL(2),NACCEL<2),VACCEL(2} 

REAL V<30>,ANET(475) 

REAL*8 DDOTSQfDGDOT 
DOOTSQ-0. 

D6D0T-V(t)*V(1)*ANET<t) 

DO 10 I>2,NQ 
NX«NXAR(I) 

7 DO 8 J«2,I 

8 DDOTSQ»DDOTSQf V < 1 ) *V ( J- 1) *ANET (NX> J- 1 ) 

1 0 D6D0TsD6D0T^V( I ) «V< I )*ANET < NX+I ) 

THE DIAGONALS ARE HANDLED SEPARATELY BECAUSE EACH OFF> 
DIAGONAL APPEARS TWICE, AND SO NUST BE DOUBLED. 
DDOTSQ-DDOTSQ^DDOTSQ+DGDOT 
DOTSQ » DDOTSQ 
RETURN 
END 
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FUNCTION DOTSQ(V,ANET) 

C CALCULATES THE INNER PRODUCT V.V RELATIVE TO THE METRIC AMET 

C 

DIMENSION LIST(30), XTMP( 8), YTMP< 7), WADJ(2A), W<26) 

1 INDEX(37), LSUBS<AC),LSUPER{39), IDADJ(23), NSYMB(36), 0U<25) 

2 PCUH(33),PRIRCM(32), CIN(22), CT0T(21), PROP(29), SPFAC<I5) 

3 0PR0P(28), V0LIN(18), V0LRT(17), DC0Nn6), PQRAT(1A>, ODEN<t?) 

4 DISS(3S), PPASS(34), PST(30), 0CINi20), PC0ND<3U ,OPRIOR(27) 

5 PAVE(13), PILE(12) 

DIMENSION VRIN(47S),8EN(999),GREF(999),ALINK(99) 

EQUIVALENCE <LINK<4I ) ,ALINK(4t ) ,GREF(8) , GEN<7) ,VRIN(7) ) 
EQUIVALENCE ( LINK(4t), LSUBS(40) ,LSUPER<39) , LIST(38), 

1 INDEX(37), NSYM6(36), DISS(35), PPASS(34), PCUM(33), 

2 PRIRCM(32), PC0ND<3U, PST(30), PR0P<29), 0PR0P(28), 

3 0PRI0R(27), U(26), 0U(2S), UADJ<24), IDADJ(23), 

4 CIN(22), CT0T(2t), 0CIN(20), 0DEN(I9), VOLIN(ld), 

5 V0LRT(17), DCON(U), SPFACdS), PQRAT(14), PAVE(13), 

6 PILEU2), XTMP< 8), YTMP( 7)) 

COMMON/CLUS/ JUNK <12), NARL , NTOP , NT6SZM ,NUANT ,L INK ( 1 4000 ) 

DIMENSION MXAR<31),LR(3),LV<3) 

EQUIVALENCE <LR< I ) ,LVRIN) , (LR(2) ,LKURT), 

1 <LR(3) ,LOVAR) , <LV< 1 ) ,iSUM) ,<LV<2) ,LSKEU) ,<LV<3) ,LOSUH) 

C 

COMMON /MISC/ HQ,MM,LR,LV,NINCLS,MXAR,UTINIT,KROOT,£PS,DELT, 

1 AMQ,0DC0N,X0VFL0,XUNFL0,UADJIN,EL1MTH,SEPTH,VFAC,AMN,SBLTH, 

2 INDXVL,UFAC,NPTSO,PORATN,SPMVTH,DUFAC,GRACTN,ANOFAC, 

3 AMOMIN ,AMOMAX , AMORAT , VOLLIM ,BIAS ,PJOIN ,VRJOIN ,US1M ,UDELSN , 

4 BETTER, MODE 

C ( END OF STANDARD BLOCK. ) 

C 

DIMENSION NTB(32) 

C 

COMMON /STPAR/UAIT,CONLV,SKBND,SKCHI,TRBND,TRCHI,URKBND,URKCHI, 

1 PACCEL<2),MACCEL(2),VACCEL<2) 

REAL V<30),AMET<475) 

00TSQ>V<1)*V<1)*ANET<1) 

DO 10 I»2,NQ 
SQaO. 

MX>MXAR(I) 

7 DO 8 Js2,I 

8 SQ»SQ+V(J-1)=»AMET(MX+J-1) 

10 DOTSQ=DOTSQ+V<I)*<<V(I)=*AMET(HX+I)+SQ)+SQ) 

C THE DIAGONALS ARE HANDLED SEPARATELY BECAUSE EACH OFF- 

C DIAGONAL APPEARS TUICE, AND SO MUST BE DOUBLED. 

RETURN 

END 


D0T00010 
DOT00020 
DOT00030 
,DOT00040 
,DOT00050 
,D0T00060 
,DOT00070 
,DOT00080 
DOT00090 
DOT00100 
DOT00110 
DOT00120 
D0T00130 
DOT00140 
D0T0015O 
DOTOOUO 
D0T001 70 
DOT00180 
DOT00190 
D0T00200 
D0TO0210 
DOT00220 
D0TO0230 
D0T00240 
D0TO0250 
DOT00260 
D0T00270 
DOT00280 
D0T00290 
D0T00300 
DOT00310 
D0T00320 
D0T00330 
D0T0034O 
D0TO0350 
DOT00360 
D0T00370 
DOT00380 
D0T00390 
BOT00400 
D0T00410 
DOT0O42O 
D0TO0430 
DOT00440 
D0T00450 
DOT00460 
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ORl^NAL PACE « 

^ QUALITY 


DOrSQA 

*» THIS ROUTINE UAS WRITTEN BY EL0GIC,INC., N0V,I?79, TO REPLACE THE DOTOOOlO 
CORRESPONDING FORTRAN ROUTINE IN THE CLASSY SYSTEH. FOR ADDITIONAL DOT00020 
INFORNATION CONTACT ELOGIC,IHC. DOT00030 

>f DOT00040 

>»FUNCTION DOTSO(V,AHET): CALCULATE THE INNER PRODUCT DOTOOOSO 



6BLC 

SP 

ADD - SUBTRACT PRECISION 

DOT00040 

2P 

SETC 

^ER' 

ER SINGLE PRECISION; DR- DOUBLE 

PREC DOT00070 

DOTSO 

CSECT 



DOTOOOSO 


USING 

DOTSO, IS 


DOT00090 


STN 

2, 12, SAVE 


DOT00100 


LH 

1,2, 0(1) 

ARCS ADDR 

DOTOOnO 


LA 

4, 4(0,0) 

I INCREHENT IS 1 

DOT00120 


EXTRN 

HISC 


D0T00130 


L 

11,>A(HISC) 


D0T00140 


USING 

HISCD,n 


D0T00150 


L 

5, HQ 

I COHPARAND IS HQ 

D0T00140 


SLA 

5,2 

HULTIPLY BY 4 

D0T00170 

!* 




DOT00180 


LR 

4,4 

J INCREHENT IS 1 

DOT00190 


LA 

8, 8(0,0) 


D0T00200 


SR 

1,8 

INDEX OF V STARTS FROH 1 

D0TO0210 


SR 

2,8 

INDEX OF AHET STARTS FROH 1 

DOTO0220 

!» 




DOT00230 


LE 

0,4(1, 4) 

V(1) 

00T00240 


HER 

0,0 

V(1)*V(1) 

O0TO02S0 


HE 

0,4(2, 4) 

♦AHET(I) 

D0TO0240 


SIP 

2,2 

0. 

DOT00270 

!» 




D0TO0280 


LR 

7,8 

DO LOOP INDEX I STARTS FROH 2 

DOT00290 

II 

LIP 

4,2 

SQ«0. 

DOT00300 


L 

12,HXAR-4(7) 

HX«HXAR(I) 

D0TO0310 


SLA 

12,2 

HULTIPLY BY 4 

D0T00320. 


AR 

12,2 

ADD ADDR AHET 

DOT00330 





DOT00340 


LR 

9,8 

DO LOOP INDEX J STARTS FROH 2 

D0TO03S0 

12 

LE 

4, 0(1,9) 

V(J-1) 

D0T00340 


HE 

4,0(12,9) 

V(J-1)<'AHET(HX+J-1) 

DOT00370 


ASP 

4,4 

SQ 

DOT00380 


BXLE 

9,4,12 


D0T00390 

nilTAA A/\A 


D0T00400 
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0«l(iJf>fAL PAGt :b 

OF POOR QUALITY 


•4 



LE 

6,U),7> 

V(I) 

DOT00410 


HE 

6,4(12,7) 

V(I)«AHET(HXH) 

DOT00420 


AlP 

6,4 

4SQ 

00T00430 


AtP 

6,4 

^SQ 

D0TO0440 


HE 

6, 4(1, 7) 

*V(I) 

D0TO0450 


AlP 

0,6 

DOTSQ 

D0TO0460 


BXLE 

7,4,11 


DOT00470 

i» 




DOTOD480 


LH 

2, 12, SAVE 


DOT004?0 


BR 

14 


DOT00500 

:* 




D0TO0510 

SAVE 

DS 

12F 


DOT00520 

4 




DOT00530 

4 




DOT00540 

4 COHNON 

/HISC/ 


D0TO0550 

NISCD 

OSECT 



D0TO0S60 

HQ 

DS 

F 


DOT00570 

NM 

DS 

F 


DOTOOS80 

LR 

DS 

3F 


D0TO0590 

LV 

DS 

3F 


DOT00600 

NINCLS 

DS 

F 


OOT00610 

NXAR 

DS 

31F 


D0TO0620 


END 



DOT00630 
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ORIGINAL PAGE IS 
OF POOP QUALITY 




c 

c 


c 


SUBROUTINE HPVS(AN,C,V) 

SETS AN»AN+V»V*C (TENSOR PRODUCT) 


^“r«H»I.R,LV,NlNCLS,«XAR,yTINlT,KROOT, EPS, belt, 


REAL AN(475>,U<30) 

L0C>0 

DO to I»1,HQ 
DO to J»t,I 
LOC^LOCft 

to AH(LOC)*AN(LOC)+U(I)*V( J)*C 
RETURN 
END 


NPVOOOtO 

HPy00020 

NPU00030 

HPV00040 

HPV00050 

NPV00060 

HPV00070 

HPVOOOBO 

HPUOOOPO 

HPMOOtOO 

HPMOOttO 

HPV0Ot2O 

NPV00130 

HPVOOtAO 

HPVOOtSO 

NPVOOtAO 

NPVOOt 70 



original pa« » 

OF POOR QUAUTf 


page 20 


SUBROUTINE MPVS(AH,C,V) NPV00010 

C SETS (TENSOR PRODUCT) NPV00020 

C NPV00030 

COHHOM /HISC/ Ha,MN,LR,LV,NINCLS,NXAfi,«TINIT,KROOT,EPS,DELT, MPV00040 

t ANQ,ODCON,XOUFLO,XUNFLO,UADJINfEUHTK,SEPTHfUFRC,RNH,SBLTH, npvoooso 

2 INOXUL,UFACFNPrsO,PORATH,SPHUTH,OUFAC,GRACTN,ANOFAC, NPU00060 

3 ANOHIN, AN0NAX,AN0RATfV0LL1N, BIAS, PJO1N,URJOINfUSIN,U0ELSN, NPV00070 

4 BETTER, NODE NPM00080 

C HPV00090 

REAL AH(475),V(30) HPV00100 

L0C«0 HPV001IO 

DO 10 1*1, HQ HPV00120 

CV»C*V(I) NPV00130 

DO 10 J»1,I HPV00140 

L0C*L0C<-1 HPV00150 

10 AHaOC)«AHaOC)fV<J)«CU HPUOOUO 

RETURN HPV00170 

END HPU00180 



J 


pags 


page 2,1 


HPVSA 

*» THIS ROUTINE MAS URITTEN BY ELOGIC,INC., N0V,I?79, TO REPLACE THE HPV00010 
« CORRESPONDING FORTRAN ROUTINE IN THE CLASSY SYSTEH. FOR ADDITIONAL NPV00020 


« INFORHATION CONTACT EL06IC.INC. 

HPV00030 





HPV00040 

>* SUBROUTINE HPVS(AN.C.V) i SETS 

AH«AH«V*V*C 

HPV00050 

NPVS 

CSECT 



HPV00040 


USING 

HPVS,15 


NPV00070 


STN 

2, 12, SAVE 


MPV00080 


LN 

1,3,0<1) 

ARGS ADDR 

HPV00090 


LA 

4, 4(0,0) 

I INCREMENT IS 1 

HPVOOlOO 


EXTRN 

NISC 


HPV00110 


L 

12,>A(HISC) 


HPV00120 


USING 

NISCD,t2 


HPVOOl 30 


L 

5, HQ 

COMPARAND IS HQ 

HPV00140 


SLA 

5,2 

MULTIPLY BY 4 

HPVOOl 50 


SR 

5,4 

SUBTRACT 4 

HPVOOlOO 

!» 




HPV00i70 


LR 

8,4 

J INCREMENT IS 1 

HPVOOl 00 


SR 

7,7 

LOG 

HPVOOl 90 

* 




HPV00200 


SR 

9,9 

DO LOOP INDEX I 

HPV00210 


LE 

2,0(2) 

C 

HPV00220 





HPV00230 

11 

LER 

6,2 

C 

HPV00240 


HE 

6, 0(3, 9) 

C«V(I) 

HPV00250 

« 




HPV002I&0 


SR 

11,11 

DO LOOP INDEX J 

MPV00270 

12 

LER 

0,6 

C«V(I) 

HPVOO20O 


HE 

0,0(3,11) 

G*V(I)*V(.') 

HPV00290 


AE 

0,0(1, 7) 

ADD AH(L£;e) 

HPV00300 


STE 

0, 0(1,7) 

STORE AH(LOC) 

HPV00310 


AR 

7,4 

ADD 1 TO LOC 

HPV00320 


BXLE 

11,8,12 


HPV00330 


BXLE 

9,4,11 


HPV00340 

III 




HPV00350 


LH 

2, 12, SAVE 


HPV00360 


BR 

14 


HPV00370 

SAVE 

DS 

12F 


HPV00380 

CONNON 

/HISC/ 


HPV00390 

NISCO 

DSECT 



HPV00400 

NQ 

DS 

F 


HPV00410 


END 



HPV00420 
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SUBROUTINE VNTV<MR,RN£T,UB) UNT00010 

C SETS UR»AH£T*VB VHT00020 

CANNON /NISC/ NQ,NHfLR,LV,NJNCLS,NXRR,UTINIT,HROOT,EPS,D£LT, VHT00030 

1 RNQ,OOCON,XOUFLO,XUNFLO,UAIIJIN,EL1NTH,SEPTH,VFAC,<IH^,SBLTH, VNT00040 

2 INDXVL,WFAC,NPTSO^PflRRTH,SPNVTH,D«FAC,GRACTt1,AMOFAC, VHT00050 

3 ANOHIN, RN0NAX,AN0RAT,V0LLIN, BIAS, PJOlN,MRJOIN,U^IN,;:c;tLSN, UNT00060 

4 BETTER, NODE, C0RLEN,SPC0R yHT00070 

REAL VA(30),VB(30),ANET(475) VHT00080 

L0CA»0 VHT00090 

00 20 I»1,HQ UNTOOtOO 

SUN-0. VNTOOttO 

00 10 J>I,I UNTO0I20 

L0CA-L0CA>l VNTO0I30 

to SUN-SUNf ANETaOCA)«UB(J) ONTOOMO 

IFU.EQ.NQ) GO TO 20 VHT00150 

JS-I+I VNT00160 

LOCB-LOCA+I VNT00170 

DO n J>JS,NQ UNT00I80 

SUN-SUN+ANET ( LOCB ) *VB ( J ) VNTOO 1 90 

tl L0CB-L0CB4J UNT00200 

20 UAm-SUN VNT00210 

RETURN VNT00220 

END UNT00230 
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V«TVA 


>» THIS 

ROUTINE UAS URI7TEN BY ELOGIC, INC., N0V,)?77, TO REPLACE 

VHT00010 

>* CORRESPONDING FORTRAN ROUTINE IN THE CLASSY SYSTEH. FOR ADDITIONAL 

VHT00020 

>» INFORHATION CONTACT ELOGIC, 

INC. 

VHT00030 

>» 




VHT00040 

subroutine ' 

VHTV(VA,AHET,VB): 

SETS VA»AHET*VB 

VHT00050 

vntv 

CSECT 



VHT00060 


USING 

VHTV,15 


VHT00070 


STH 

2, 12, SAVE 


VHT00080 


LN 

t,3,0(t) 

ARGS ADDR 

VHT00090 


LA 

4,4<0,0> 

I INCREHENT IS 1 

VHT00100 


EXTRN HISC 


VHTOOItO 


L 

7,>A(HISC) 


VHT00120 


USING 

NISCD,7 


VHT00130 


L 

5, HO 

COHPARANO IS HQ 

VMT00140 


SLA 

5,2 

HULTIPLY BY 4 

VHT00150 


DROP 

7 


VHT00160 

l» 




VHT00170 


LR 

8,4 

J INCREHENT IS 1 

VHT00180 


SR 


INDEX OF VA STARTS FROH 1 

VHT00190 


SR 


INDEX OF AHET STARTS FROH 1 

VHT00200 


SR 

3,4 

INDEX OF VB STARTS FROH 1 

VHT00210 


SR 

7,7 

INITIALIZE LOCA 

VHT00220 





VHT00230 


LR 


DO LOOP INDEX I STARTS FROH 1 

VMT00240 

It 

LE 

i,»E^O.O" 

INITIALIZE SUH 

VHT00250 

l» 




VHT00260 


LR 

11r4 

DO LOOP INDEX J STARTS FROH 1 

VHT00270 

12 

AR 

7,4 

ADD 1 TO LOCA 

VHT00280 


LE 

0,0(2, 7) 

AHET (LOCA) 

VHT00290 


HE 

0,0(3,11) 

HULTIPLY BY VB(J) 

VHT00300 


AER 

6,0 

ADD TO SUH 

VHT00310 


BXLE 

11,8,12 


VHT00320 

!» 




VMT00330 


CR 

9,5 

CON. ^RE INDEX 1 UITH ITS COHPARANO 

VMT00340 


BZ 

A12 


VHT00350 


LR 

n,9 

I 

VHT00360 


AR 

11,4 

JS«If1 

VHT00370 


LR 

12,7 

LOCA 

VHT00380 


AR 

12,9 

LOCA *I 

VHT00390 

13 

LE 

0,0(2,12) 

AHET(LOCB) 

VHT00400 


HE 

0,0(3,11) 

HULTIPLY BY VB(J) 

VHT00410 


AER 

‘6,0 

ADD TO SUH 

VHT00420 


AR 

12,11 

ADO J TO LOCB 

VHT00430 


BXLE 

11,4,13 


VHT00440 

* 




VHT00450 

A12 

StE 

6, 0(1, 9) 

VA(I) 

VHT00460 


BXLE 

9,4,11 


VHT00470 

!* 




VHT00480 


LN 

2, 12, SAVE 


VHT0049C 


BR 

14 


VHT00500 

SAVE 

DS 

12F 


VHT00510 

■* COHHON 

/HISC/ 


VHT0052C 

NISCD 

DSECT 



VHT00530 

ilfl 

DS 

F 


VHT00540 


END 



VHT00550 
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SUBRQUTIHE VP\;(VA,FAC,gB) VPV00010 

C SETS VA»VA+FAC*VB Vpy00020 

c vpyooo3o 

COMMON /MISC/ MQ,NM,U»l-V,NIMCLS,MXAR,UTINrr,KROOT,EPS,DELT, ypyoOOAO 

1 AMO,ODCON,XOVFLO,XUNFLO,HADJIN,ELIMTH,SEPTH,VFAC,AMM,SBLTH, VPM00050 

2 INDXVL,UFAC,NPTSO,PQRATH,SPMVTH,DUFAC,GRACTM,AHOFAC, VP«JOOO^O 

3 AMOM1N,AMOMAX,AMORAT,VOLLIM,BIAS,PJOIN,VRJ01N,USIM,UDELSH, Vpg00070 

4 BETTER, MODE, C0RLEN,SPC0R VPV00080 

c vpyooopo 

REAL VA(30),VB(30) VPyOOtOO 

00 10 I»1,M0 VPMOOilO 

10 VA<I)*VAm+VBa)*FAC Vpy00120 

RETURN VPM00130 

END VPMOOHO 


ORIGINAL PAGE IS 
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=» THIS 

ROUTINE UAS URITTEN BY 

ELOGIC, INC., NOV, 1979, TO REPLACE THE 

VPV00010 

CORRESPONDING FORTRAN ROUTINE IN THE CLASSY SYSTEH. FOR ADDITIONAL 

VPV00020 

'» INFORNATION CONTACT ELOGIC, 

INC. 

VPV00030 

•* 




VPV00040 

SUBROUTINE ' 

VPV<VA,FAC,VB): 

VA(I)«VA(I)+FAC*VB(I) 

VPV00050 

VPV 

CSECT 



VPV00060 


USING 

VPV,15 


VPV00070 


STN 

2, 7, SAVE 


VPV00080 


LN 

t,3,o<n 

ARGS AODR 

VPV00090 


LE 

2,0(2) 

ARG2 

VPV00100 


LA 

4, 4(0,0) 

INCREHENT IS 1 

VPV00110 


EXTRN 

HISC 


VPV00120 


L 

7,«A(NISC) 


VPU00130 


USING 

HISCD,? 


VPV00140 


L 

S,HQ 

COHPARAND IS HQ 

VPV00150 


SLA 

5,2 

HULTIPLY BY 4 

VPV00160 


SR 

5,4 

SUBTRACT 4 

VPV00170 

.•» 




VPV00180 


SR 

2,2 

DO LOOP INDEX I STARTS FROM 1 

VPV00190 

IC 

LER 

0,2 

FAC 

VPV00200 


HE 

0,0(3, 2) 

FAC*VB(I) 

VPV00210 


AE 

0,0(t,2) 

VA(I) 

VPV00220 


STE 

0,0(1, 2) 

STORE VAU) 

VPV00230 


BXLE 

2,4, IC 

BRANCH LOCATION 

VPV00240 





VPV00250 

• 

LN 

2, 7, SAVE 


VPV00260 


BR 

14 


VPV00270 

SAVE 

DS 

7F 


VPV002BO 

« CONHON 

/NISC/ 


VPV00290 

HISCD 

DSECT 



VPV00300 

HQ 

DS 

F 


VPV00310 


END 



VPV00320 
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Appendix 2: Program listings for the fast exponential system. 

This appendix contains the routines used by the fast expon- 
ential system. They include the routine XP, the table generator 
XPKEP, the modified version of CLINIT, showing where XPREP is 
called; and the modified version of STATIS, showing the calls to 
XP. (The version of STATIS displayed also contains modifications 
for the additional statistics collection, which will be the 
siobject of an additional report.) The modification to CLINIT 
consists only, of the one line calling XPREP, which may be exe- 
cuted even if XP is not used. Elogic recommends that the version 
containing the call to XPREP be used at all times. 
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THIS ROUTINE UAS WRITTEN BY ELOGIC.INC., NOV, 1 979, TO REPLACE THE 

XP 

00010 

CORRESPONDING FORTRAN ROUTINE 

IN THE CLASSY SYSTEH. FOR ADDITIONAL 

XP 

00020 

it- INFORNATION CONTACT ELOGIC,INC. 

XP 

00030 




XP 

00040 

EXPONENTIAL 

ROUTINE, EXP(-X/2) 


XP 

00050 

XP CSECT 



XP 

00060 

USING 

*,15 


XP 

00070 

L 

1,0(1) 

ARG ADDR 

XP 

oooso 

LE 

0,0(1) 

ARC 

XP 

00090 

AD 

0,SHIFURD 

EXTRACT HANTISSA 

XP 

00100 

STD 

0, CARVE 


XP 

00110 

L 

1, CARVE 

GET UPPERUORO 

XP 

00120 

S 

1,UPCARVE 

CHK ARG BETWEEN 

XP 

00130 

BNZ 

NSEXP 

0 AND 16 

XP 

00140 

XPSEQ IC 

l,CARVE+4 

GET UPPER BYTE 

XP 

00150 

HVI 

CARVE+4,X"3F" 

PUT IN EXPONENT 

XP 

00160 

LE 

0,CARVEH 

LOAD XI 

XP 

00170 

AE 

0,XPD 

CALC EXPRESSION 

XP 

00130 

ALR 

Ifl 

INDEX 

XP 

00190 

UR 

1,1 


XP 

00200 

HE 

0,XPAR(1) 

EXPONENTIAL TABLE 

XP 

00210 

BR 

14 

(BYE) 

XP 

00220’ 

i» 



XP 

00230 

NSEXP C 

1,H2NCARV 

CHK UNDERFLOW 

XP 

00240 

BH 

ZEROIT 


XP 

00250 

C 

1,HAXCARV 

CHK OVERFLOW 

XP 

00260 

BP 

OUTRANGE 


XP 

00270 

SR 

1,1 


XP 

00290 

IC 

1,CARVE^4 

(DUP XPSEQ) 

XP 

00290 

HVI 

CARVE+4,X''3F" 


XP 

00300 

LE 

0,CARVE+4 


XP 

00310 

AE 

0,XPD 


XP 

00320 

ALR 

1,1 


XP 

00330 

ALR 

1,1 


XP 

00340 

HE 

0,XPAR(1) 


XP 

00350 

L 

1, CARVE 

GET BIG PART 

XP 

00360 

ALR 

1,1 

INDEX -USES SPEC 

XP 

00370 

ALR 

1,1 

PROPERTIES OF SHIFURD 

XP 

00330 

HE 

0,XPBIG(1) 

BIG EXPONENTIAL TABLE 

XP 

00390 

BR 

14 

(BYE) 

XP 

00400 




XP 

00410 

ZEROIT LE 

0,=F-'0-' 

UNDERFLOW - 

XP 

00420 

BR 

14 

RETURN 0. 

XP 

00430 

* 



XP 

00440 

OUTRANGE L 

1,0UTCNT 


XP 

00450 

A 

1,=FU'' 


XP 

00460 

ST 

1,0UTCNT 


XP 

00470 

B 

ZEROIT 


Al' 

00430 
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t» 




XP 

00490 

If 




XP 

00500 

CARVE 

DS 

D 


XP 

00510 

SHIFURD 

DC 

A(XM78000000 

SHIFTS ARGUHENT TO GET 

XP 

00520 


DC 

F'O' 

HANTISSA 

XP 

00530 

UPCARVE 

EQU 

SHIFURD 

ELIHS SHIFURD 

XP 

00540 

HINCARV 

DC 



XP 

00550 

HAXCARV 

DC 

?' 2 \' 


XP 

00560 

If 




XP 

00570 

XPC0« 

DS 

OD 


XP 

00580 


ENTRY 

XPCOH 


XP 

00590 

OUTCNT 

DC 

F"0" 


XP 

00600 

XPXTRA 

DC 

F-'O^ 


XP 

00610 

XPD 

DS 

D 


XP 

00620 

XPAR 

DS 

256F 


XP 

00630 


DS 

21F 

XPDIG NEG ARG 

XP 

00640 

XPBIG 

DS 

22F 


XP 

00650 


END 



XP 

00660 
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SUBROUTINE XPREP XPROOOlO 

COHHON /XPCOH/ 0UTCNT,XPXTRA,XPD,XPAR0(256),XPBIG0(43) XPR00020 

REALMS XPD, ALPHA, C2,Dt XPR00030 

XPR00040 

BETA* 1./ 16. XPR00050 

ALPHA-BETA/4. XPR00060 

C2»2.*ALPHA XPR00070 

D1»3.*(DSINH<ALPHA)-ALPHA*BC0SH(ALPHA))/(C2*ALPHA*ALPHA> XPR00080 

XPD*DSINH(ALPHA)/<ALPHA*D1) - C2 XPR00090 

DO 10 N*1,2S6 XPRD0100 

10 XPAR0(N)=D1*DEXP<-C2*(N-1)-ALPHA) XPR0O1I0 

DO 20 1*1,43 XPR00120 

20 XPBIG0(I)*DEXP(-8.D0*a-22)) XPR00130 

RETURN XPR00140 

END XPR00150 
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POOR 


fS 

QiJALl-a 


SUBROUTINE CLINIT(KROT) 

THIS ROUTINE CONTAINS THE VARIOUS STATEMENTS NECESSARY TO 
INITIALIZE THE CLUSTERING ALGORITHM. 

DIMENSION LIST(38), XTMP( 8), YTMP< 7 ), UADJ(24), W(26) 

t IN0EX(37), LSUBS<40),LSUPER(3?), IDADJ<23), NSYHB(36), 0U(2S) 

2 PCUM(33),PRIRCM(32), CIN(22), CT0T(2I), PR0P(29), 3PFAC(15) 

3 0PR0P(28), VOLINdS), V0LRT(I7), DC0N(t6), PQRAT(H), ODENd?) 

4 DISS<35), PPASS(34), PST(30), 0CIN(20), PCQND(31) ,OPRIQR (27) 

2 PAVE(13), PILE(12) 

DIMENSION VRIN(475),GEN(999),GREF(999),ALINK(99> 

EQUIVALENCE (LINK(4t ) ,ALINK(41) ,GREF(8) , GEN(7) ,VRIN(7) ) 
EQUIVALENCE ( LINK(41}, LSUBS(40),LSUPER(39) , LIST(38), 

1 INDEX(37), NSYMB(36), DISS(35), PPASS(34), PCUM(33), 

2 PRIRCM(32), PC0ND(3t), PST(30), PR0P(29), 0PR0P(23), 

3 0PRI0R(27), U(26), 0U(2S), UADJ(24). IDADJ(23), 

4 CIN<22), CT0T(21), 0CIN(20), 0DEN(I9), VOLINdS), 

5 V0LRT(17), DC0Nd6), SPFACdS), PQRAT(14), PAVEd3), 

6 PILE(12), XTMP( 8), YTMP( 7)) 

COMMON/CLUS/ JUNKd 2) ,NARL,NTOP,NTBSZM,NUANT,LINK(I4000) 

DIMENSION MXAR(31),LR(3),LV(3) 

EQUIVALENCE (LRd ) ,LVRIN) , (LR<2) ^LKURT) , 

1 <LR(3),L0VAR),(LVd),LSUM),(LV(2),LSKEU),(LV(3),L0SUM) 

C 

COMMON /MISC/ MQ,MM,LR,LV,NINCLS,MXAR,UTINIT,KROOT,EPS,DELT, 
t AMQ,ODCON,XOVFLO,XUNFLO,UADJ1N,ELIHTH,SEPTH,VFAC,AMN,SBLTH, 

2 INDXVL,UFAC,NPTSO,PQRATH,SPMVTH,DUFAC,GRACTM,AMOFAC, 

3 AM0MIN,AM0MAX,AM0RAT,V0LLIM,6IAS,PJ0IN,VRJ0IN,USIM,UDELSM, 

4 BETTER, MODE 

C ( END OF STANDARD BLOCK. ) 

C 

DOUBLE PRECISION XTEHP,YTEMP,ZTEHP,DURK,DURKD 
C 

COMMON /STPAR/UAIT,CONLV,SKBND,SKCHI,TRBND,TRCHI,URKBND,URKCHI, 

1 PACCEL(2),MACCEL(2),VACCEL(2) 

C 

COMMON/CLUSTR/ IBEGIN,TOTURD,CLSNAM,IPT,NOFLD, SYM(61) , 

1 LNCAT, PRNT(4) , KLBC , PRIME, PROUT, TOTPIX, 

2 SCRAM1 ,BUFPIX,BUFTOT,MBUFSD,NDUMP,LBUFD 

3, MAXBF, AREA, NUDS, NUDRS, NPTS, LBUF, IQ1,N0CYCL, NCL 
C 

INTEGER TOTURD,SYM,PRNT,PRTME, PROUT, TOTPIX, SCRAM1,BUFPIX,BUFT0T 
1 ,CLSNAM 
C 

COMMON /MXLL/ MXLLUT, MXLLFN, RELPRP(200) 

C 

COMMON /INITL/UTNEU,DEVINI,CHANIN 
CHIVAL(DF)=DF*d.-.222/DF+C0NLV*SQRT(.222/DF))*=»3 
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CLI00420 
CLI00430 
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C 

C 

C 


C 

C 


INITIALIZS ADDRESS 
CALL DATFIX 
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OP POOR QUALrry 


AHQ»HQ 

HQS>HQ«HQ 

DEFINE VALUE OF SEPTH IN TERNS OF CHI SQUARE VALUE 
DFT = ANQ + 1 
SEPTH = (CHIVAL(DFT))/2 


C 


UE FIRST SET UP VARIOUS INDEX ARRAYS FOR A PARTICULAR 
NUHBER OF CHANNELS NQ. 

SET UP THE TRIANGULAR POSITION ARRAY HXAR. 

NN=0 

DO 10 1=1,31 
NXAR(I)=NH 




10 NH=HH+I 

NH=NXAR(HQ^1) 

AHN=NN 

SET UP TABLES FOR THE XP FUNCTION 
CALL XPREP 

NON UE SET UP THE ORIGIN VECTORS, LR AND LV, OF THE VARIOUS ARRAYS 
AND VECTORS IN A CLUSTER NODE. 

NINCLS=1 

THIS CONSTANT NUST BE SET TO THE NUHBER OF ARRAYS 
DO 21 1=1,3 
LR(I)=NINCLS 

21 NINCLS=NINCLS+HH 
DO 22 1=1,3 
Ly(I)=NINCLS 

22 NIHCLS‘NINCLS+HG 


NSCALS = 3S 

NINCLS=NINCLStNSCALS-l 

UE NUST ALSO SET UP SOME THRESHOLDS FOR USE BY THE STATISTICAL 
SYSTEH. 

SKCHI = <AHQ+2.)*(AHQ-^4.)«CHIVAL<AHQ) 


URKCHI=ANQ*(AHQ+4.)*(AHQ+6.)/(AHQ-.99?)*CHIVAL<AHH-1.) 

TRCHI=C0NLV*C0NLV*<AH0*<AHQ+2.)*(AHQ+3.)*8.) 

UE NOU CREATE THE HEAD NODE OF THE CLUSTER TREE. THIS IS NOT 
AN ACTUAL CLUSTER, AND DOES NOT HAVE STORAGE FOR ANY 
OF THE STATISTICAL ARRAYS. 

NPTS0=0 

KROT=HORSTR<NSCALS) 

HAKE FIRST NODE START AT AN ODD NUHBER 
IF (H0D<NT0P,2) .NE. 1 ) • NTOP * NTOP + 1 
LINK(KR0T)=-262139 
LSUPER(KR0T)=-262I42 
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CLI00620 
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CLI00750 

CLI00740 
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CLI00850 
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c 


c 


c 


IDADJ(KR0T)»999999 

INDEX(KR0T)=-0 

SPFAC(KR0T)»99999. 

M(KROT)=UTINIT 

OU<KROT)=y(KROT) 

PQRAT(KR0T)=0. 

PROP(KROT)*!. 

0PR0P(KR0T)»l . 

CIN(KR0T)»y(KR0T) 

0CIH(KR0T)=CIN(KR0T) 

CT0T(KR0T)=0. 

0DEN(KR0T)=U(KR0T) 

PRIRCH i KROT ) - 1 • 

NEXT THE INITIAL NODE IS SET UP, TOGETHER WITH SOHE CONTROL 
ODCON -60. 

PROPI = 1. 

INDXVL = 0 

57 KFIR = NEUCLS(KROT, PROPI, WTINIT) 

LSUBS(KROT) -KFIR 
LINK(KFIR)=0 

DEV2UT >DEVINI«DEVINI«WTINIT 
LA=NORSTR(HQS} 


CLI00970 
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CLI00990 
CLI01000 
CLI01010 
CLI01020 
CLI01030 
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CLI01070 
CLI0I080 
CLI01090 
THRESHOLDCLI01IOO 
CLI01110 
CLI01120 
CLI0H30 
CLIO1M0 
CLI01150 
CLI01I60 
CLI01170 
CLIO1180 


DO 52 I»1,NQ 

GREF(KFIR+LSUH+I)=UTINIT#CHANIN 
DO 52 J>1,NQ 
IJ=J-1+(I-1)*HQ+LA 
ALINK(IJ) - 0 . 

IF(I.EQ.J) ALINK(IJ)=>DEV2UT 
52 CONTINUE 

CALL ZER0CL(KFIR,ALINK(LA),GREF(KFIR+LSUN+1),6REF(KFIR+LSKEU+1), 

1 GREF(KFIR+LKURT+1), GREF (KFIR+LOSUH+t ) ,GREF (KFIR+LOVAR+1 ) ) 

C/;LL FREE < LA, NOS) 

SET SWITCHES FOR NAX LIKLIHOOD LABELING 
NXLLWT = 0 
UADJ(KFIR) >UADJIN 
NXLLFN = 23 

PRINT 273,HQ,C0NLV,TRCHI,SKCHI,URKCHI 
273 FORMAT ('1 CONFIDENCE LEVELS', 14,^ CHANNELS', FB. 4,' CHISQUARES', 
1 3F12.2) 


INITIALIZE ADDRESS : 

CALL STATAD(6ENaSUH),GEN(LSKEU),GEN(LKURT), 
1 6£N(L0SUH),GEN{L0VAR)) 

RETURN 

END 
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